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^ ■ ABSTRACT 

\ We study the redshift space correlation function of galaxy clusters for observational 

^ ■ samples constructed in different surveys. We explore correlation amplitudes, pairwise 

' velocity distributions and bias factors. Systematics in cluster identification procedures 

T-H , are the main source of biased estimates of the correlation amplitude and inferred ve- 

CO ' locity dispersions. We find that the large elongations along the line of sight in the Abell 

O^l I catalogue cannot be explained solely in terms of the errors in distance measurement 

, originating from using a small number of galaxies. The inclusion of a significant frac- 

■ tion of galaxies and systems not physically bound to the clusters are responsible for 

\ this large systematic effect. We also find a significant dependence of the redshift space 

f~i i distortion of the correlation function on the cluster Bautz-Morgan type, an effect that 

C^h' may rely on the fact that, due to the regular appearance of low Bautz-Morgan type 

Q I clusters, a sample of such objects would be less contaminated. We confirm that the 

5^ ■ effect of a low number of redshift measurements, riz , is to increase the redshift-space 

' correlation length and bias factor. The results are very stable for > 10. We also test 

. the effects of different in the catalogue of groups/clusters derived from the Updated 

IL" ' Zwicky Catalogue and find that, for cluster samples identified in redshift surveys, even 

5^ \ a low number of redshifts ~ 5, is sufficient to provide reliable results. By compar- 

^ . ing our results with those of numerical simulations we explore the strong influence 

' on the clustering distortion pattern in redshift space from effects associated with the 

cluster identification procedure from two dimensional surveys. The identification of 
clusters in X-ray surveys improves this situation, although there are still systematic 
effects which are probably due to identification of optical sources in the determination 
of cluster redshifts. These systematics are particularly strong for the most luminous 
X-ray selected clusters in the Extended Bright Cluster Survey, which exhibits very 
large anisotropies, comparable to those present in the Abell catalogue. Our results 
demonstrate that forthcoming large redshift surveys will be extremely important for 
the construction of new samples of groups and clusters as well as improving the deter- 
mination of optical and X-ray selected cluster distances, essential for reliable analyses 
of the large scale structure. 

Key words: methods: statistical - methods: numerical - large-scale structure of 
Universe - galaxies: clusters: general 



1 INTRODUCTION 

Clusters of galaxies can be used as tracers of the large scale 
structure of the Universe. However, their identification is 
not free from systematics, for example the Abell catalogue 



is known to be afflicted by projection effects (Sutherland 
1988; Sutherland & Efstathiou 1991; see also Lucey 1983), 
which are somewhat reduced when restricting the sample of 
clusters to high richness class (Milier et al. 1999). 



© 0000 RAS 



2 A^. D. Padilla & D. G. Lambas 



Cluster peculiar motions produce an apparent distor- 
tion of the clustering pattern as measured by the two-point 
correlation function in redshift space, {(a, tt) (Croft & Efs- 
tathiou 1994; Bahcall, Cen & Gramann 1994), with tr and 
TT being the separations perpendicular and parallel to the 
line of sight, respectively. The non-lineax virialized regions 
induce elongations along the line of sight which make it pos- 
sible to estimate the one-dimensional pairwise rms velocity 
dispersion, < >^'^ (Davis & Peebles 1983). The infall 
onto ovcrdcnso regions in the form of bulk motions dom- 
inates at large scales resulting in a compression of the ^ 
contours along the direction of the line of sight. 

The presence of spurious large anisotropics along the 
line of sight of Abell clusters has been largely discussed 
(e.g., Sutherland 1988, Sutherland & Efstathiou 1991, Bah- 
call & Soneira 1983; Postman, Huchra & Geller 1992). Aris- 
ing mainly from inhomogeneities in the detection of clusters 
in two dimensions (ic. prone to strong projection effects), 
these observed large anisotropics arc extremely unlikely to 
exist for systems of galaxies within a hierarchical scenario 
of structure formation. Therefore, the clustering signal along 
the line of sight in the Abell catalogue could be artificially 
enhanced (Sutherland 1988; Sutherland & Efstathiou 1991; 
see also Lucey 1983) due to significant projection effects 
which could bias the determination of the correlation func- 
tion as well as the mass function of these systems. A velocity 
broadening of 2000kms~'^, obtained by Bahcall, Soneira & 
Burgctt (1986) from their analysis of the clustering of Abell 
clusters, was interpreted as arising from a combination of 
cluster peculiar motions as well as geometrical distortions 
of superclusters. However, van Haarlem ct al. (1997) find 
that a third of Abell clusters of richness class _R > 1 would 
not be real physically bound systems but simply projections 
of galax;ies and groups along the line of sight. It has to be re- 
called that early surveys of Abell clusters contain large frac- 
tions of low richness clusters (Abell richness class _R = 0), 
which were not intended to form complete samples suitable 
for statistical analyses. Miller et al. (1999) analysed samples 
of Abell clusters with richness R > 1, with new accurate 
determinations of cluster positions using several galaxy red- 
shifts. These authors find that the clustering signal along 
the line of sight is greatly reduced when compared with the 
Bahcall & Soneira (1983) results. The anisotropy is further 
reduced after the orientation of two superclusters that are 
elongated along the line of sight is changed. More recently. 
Miller et al. (2002) found that the inclusion of clusters in su- 
perclusters has a significant incidence in the anisotropies of 
the redshift-space correlation function. A possible interpre- 
tation of this result is that cluster identification in overdense 
regions is likely to be more affected by projection effects. 
Peacock & West (1992) also found that restricting attention 
to higher richness Abell clusters removed the strong radial 
anisotropy seen in the clustering measured from earlier sur- 
veys. 

A more recent generation of cluster catalogues presum- 
ably less affected by identification biases (APM: Dalton et 
al. 1992, 1994, 1997; Cosmos: Lumsden et al. 1992) was 
drawn from machine-scanned survey plates with better cali- 
brated photometry. With the aim of reducing biases in clus- 
ter richness determinations arising from projection effects, 
the search radius used to define clusters in these machine 
based catalogues is smaller than that used by Abell. As ex- 



pected, the clustering pattern found in these more recent 
cluster redshift surveys does not exhibit very large enhance- 
ments along the line of sight; furthermore, the trend of in- 
creasing correlation amplitude with decreasing space den- 
sity of clusters is weaker than that found for Abell clusters 
(Croft et al. 1997) and is similar to the trend expected in 
current models, where there is only a weak dependence of 
correlation length on cluster space density, as is obtained in 
other all sky X-ray surveys (eg REFLEX sample, Collins et 
al. 2000), and samples of Abell clusters confirmed in X-rays 
(X-ray Bright Abell Cluster Sample, XBACS, Ebehng et al. 
1996). For these samples, systematic projection effects along 
the line of sight are expected to be significantly lower than 
in the Abell catalogue (Abadi, Lambas & Muriel 1998; Bor- 
gani et al. 1999; Collins et al. 2000). However, it is important 
to consider that the analysis of high richness Abell clusters, 
performed by Miller et al. (1999), gives results comparable 
to the amount of distortion of the clustering pattern found 
for APM clusters (see Fig 5 of Miller et al. 1999). 

In this paper, we use the results of an analysis of the red- 
shift space clustering of massive dark matter haloes in mock 
cluster samples (Padilla & Lambas, 2003, Paper I hereafter) 
extracted from the tCDM Hubble Volume simulation. The 
results from Paper I are used in the assessment of projection 
biases on catalogues of clusters identified in two dimensions 
and three dimensions. In paper I wc studied statistically the 
redshift-space correlation function as a function of coordi- 
nates parallel and perpendicular to the line of sight, through 
the pairwise velocity dispersions, redshift-space correlation 
length and bias parameters obtained from difi'erent mock 
cluster samples constructed using several cluster identifica- 
tion algorithms, and different number of redshifts and clus- 
ter search radii. Spurious distortions in the redshift space 
cluster correlation function could originate in various sys- 
tcmatics associated to these issues. For instance in several 
catalogues, many cluster distances are determined by a sin- 
gle galaxy redshift, usually the brightest cluster member and 
therefore are subject to projection contamination from back- 
ground and foreground galaxies. 

In this paper we present an analysis of the correlation 
function of observational cluster data, including optical and 
X-ray cluster samples. In section 2 we present the statistical 
procedures adopted and we summarise the results obtained 
from the study of different mock samples in Paper I that 
concern the present study. These mock samples are used 
for assessing the degree to which different problems could 
contribute to the anisotropies detected in observational cor- 
relation functions in Section 3. In this section we study the 
effects of the number of redshifts used in the determination 
of cluster distances, the effects arising from selecting sam- 
ples with different Bautz-Morgan cluster types (Bautz & 
Morgan, 1970) and the different behaviour of samples with 
high and low x-ray luminosities. In section 4, we analyse 
groups of galaxies obtained from the Updated Zwicky Cat- 
alogue, hereafter UZC , which, in contrast to other sets of 
clusters studied in this paper, were identified from redshift 
space data. In section 5 we discuss our results and present 
the main conclusions. 
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2 STATISTICAL PROCEDURES AND 

RESULTS FROM MOCK CATALOGUES 

We followed the procedures described in paper I to mea- 
sure the redshift-space correlation function as a function of 
the pair separation parallel and perpendicular to the line 
of sight, a and n respectively. The iso-corrclation contours 
(curves of equal amplitude of ^(cr, tt) = ^^*^) provide a suit- 
able measurement of the pairwise velocities of galaxies by 
comparing the measured and predicted contour levels. These 
contour levels are approximated by the functions r"* (9) and 
r''{9) for the measured and predicted correlations respec- 
tively, where 9 is the polar angle measured from the direction 
perpendicular to the line of sight (See Paper I for details). 

In order to make a proper comparison between model 
predictions and observations, we adopt the real-space cor- 
relation function of the mass distribution in a ACDM cos- 
mology. We use a mass density parameter Qmatter = 0.3, 
vacuum energy density = 0.7, and a CDM shape param- 
eter r = 0.2 consistent with recent CMB results (Nottor- 
field et al. 2002). The normalisation of the power spectrum, 
178 = 0.9, is taken from recent constraints obtained from the 
2dF Galaxy Redshift Survey (Norberg et al. 2002). Given the 
relatively large uncertainties in the cluster correlation func- 
tions the adopted transfer function docs not include baryons, 
however our choice of F implicitly allows for a combination 
of baryon fraction and an effective value of F as discussed 
in Eisenstein & Hu (1996). We obtain the real-space cor- 
relation function by Fourier transforming the CDM power 
spectrum 
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and then use this to evaluate the theoretical prediction for 
^^(cr, tt). We search for the optimum values of the scale inde- 
pendent bias parameter b and W12 by minimising the quan- 
tity x', 
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where we have chosen to compare a set of discrete levels, 
/ = 0.6,0.8,1.0,1.2 and 1.4, of the redshift-space correla- 
tion function ^(cr, tt) amplitude instead of comparing values 
of the correlation function on a grid of a and tt distances. 
Our decision is based on the fact that more reliable and sta- 
ble results arc obtained when using the isoplcth comparison 
technique (Padilla et al. 2001). We notice that the effect of 
adopting a different normalisation for the theoretical power 
spectrum mainly affects the bias parameter b, leaving the 
pairwise velocity W12 mostly unchanged. 

We briefly summarise the results obtained for the sta- 
tistical properties of cluster samples drawn from mock cat- 
alogues presented in Paper I. These mock catalogues have 
been obtained from one of the mock APM galaxy catalogues 
extracted from the rCDM Hubble Volume simulation by the 
Durham Extragalactic Group (Evrard et al., 2002) following 
a procedure similar to that used in Cole et al. (1998). 

In Paper I, we used three main algorithms to search 
for groups in the rmmcrical simulations, and constructed 
four samples of clusters from the simulations. The aim was 
to understand the effects of cluster selection procedures on 
the cluster correlation functions, and in particular, to es- 



tablish the origin of the large elongations along the line of 
sight present in the measured correlation functions. Sample 
1, constructed using a friends-of- friends algorithm applied 
to the particles in the simulation box, has a lack of a sig- 
nificant number of spurious clusters. The measurements of 
anisotropies in the correlation function clearly show the ex- 
pected infall pattern. Sample 2 allows to understand the 
results obtained from the identification of groups in a mock 
catalogue with redshift information following the lines de- 
scribed by Huchra & GcUer (1982). In this case the results 
from Paper I indicate that the infall pattern is not as readily 
visible as in the results from Sample 1, but a small ~ 5 
is enough to produce it. Sample 3 was constructed following 
the prescription described by Lumsdcn et al. (1992), using 
an angular search radius rc = 1.0h~^Mpc. This procedure 
was applied to the same mock catalogue as the second al- 
gorithm, and the resulting sample was found to show severe 
projection effects. The results from this sample will be used 
to assess the degree of these effects in observational samples 
of cluster of galaxies identified from angular data. In Paper I 
we also extracted another sample combining the mock cata- 
logues defined above, in order to mimic as closely as possible 
the methods used in obtaining the observational samples. 
From the clusters identified using angular data (sample 3), 
we drawn a restricted set of clusters of galaxies, whose an- 
gular positions lie within the cluster search radius from a 
3D identified cluster, sample 4. 

As a summary of the results of our mock samples, we 
recall figures 11,12, and 13, of Paper I, where the behaviour 
of the correlation contours is assessed by measuring the pair- 
wise velocity dispersion W12 and the correlation length so as 
a function of nz- We found that in the case of the identi- 
fication of clusters from angular data, the values of so and 
W12 are higher as the underlying value by as much as an 
order of magnitude. In the case of clusters identified from 
3-dimensional information (sample 2), these quantities con- 
verge rapidly to their actual values, which are obtained when 
a small number of galaxies, not much larger than 5, is used 
to determine the cluster distance. 

For more detailed descriptions of the characteristics of 
each sample, analysing the differences in elongations arising 
from using different riz, and search radius rc, we refer the 
reader to Paper I. 



3 REDSHIFT SPACE DISTORTIONS OF 

GROUPS AND CLUSTERS OF GALAXIES 

In this section we compute the two-point correlation func- 
tion in redshift space using different cluster catalogues. 

Following usual procedures we estimate ^"(cr, tt) by 
counting pairs in the data and in a random catalogue with 
the same angular limits and radial selection function than 
the observational sample. Data-data and random-random 
pairs {Ndd and Nrr respectively) are binned as a function of 
separation in the two variables a and tt. We adopt the Davis 
& Peebles (1983) estimator, 



r(^,7r) = 



- 1 



(3) 



where no and nR are the number of data and random points 
respectively. 
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Our analyses are centred on large correlation amplitude 
signals which are not strongly sensitive to the choice of es- 
timator (see Hamilton 1993). 

In Figure 1, we show the contours levels of ^°((t, tt) 
for UZC groups (loft panel) . There are no signals of distor- 
tions along the line of sight direction for the groups consis- 
tent with low pairwise velocities of groups. By comparison 
with the results for samples of clusters identified in two- 
dimensional surveys shown in the other panels of this figure 
it is evident the existence of either large pairwise velocities 
or large projection biases in these samples, as also discussed 
by Sutherland (1988). Such elongations could originate in 
the systematic presence of groups along the line of sight, in 
the fields of clusters identified from angular data. From the 
theoretical point of view, strong elongations along the line of 
sight axe not expected in a hierarchical scenario of structure 
formation. 

This section is divided as follows: In section 3.1 we study 
the correlation function of APM and Abell clusters. We di- 
vide the Abell samples according to the number of redshifts 
used in the determination of the cluster distances in section 

3.1.1, and also according to Bautz-Morgan type in section 

3.1.2. Cluster samples identified in X-rays are studied in 
section 3.2; the results from subdividing the EBCS sample 
using a x-ray luminosity cut are shown in section 3.2.1; and 
the outcomes of sub-dividing the XBACS sample accord- 
ing to different Bautz-Morgan type are explained in section 
3.2.2. Finally, the inferred values of relative velocities and 
redshift-space correlation lengths are presented in section 
3.3, and the bias factors in section 3.4. 

3.1 Clusters samples selected from angular 
surveys 

In this section we study the iso-correlation contours ob- 
tained from cluster catalogues identified using surveys of 
gala^xies containing only angular positions. As has already 
been pointed out, studies of such samples have shown large 
elongations along the line of sight even when, in light of the 
previous discussion, we would have expected infall patterns 
from samples of massive clusters (Padilla & Baugh, 2002) . 

We measured the correlation function of APM and 
Abell clusters in order to compare them to the results from 
the UZC groups, shown in the left panel of figure 1. We 
also show in this figure the results for APM clusters (Dalton 
et al., 1992; middle panel), and Abell clusters with mea- 
sured redshifts (Struble & Rood 1999, right panel). The 
sample of Abell clusters consists only of clusters with de- 
termined Bautz-Morgan type. It can be seen that there is 
a marked difi'erence in the distortion pattern between the 
UZC groups, which show an infall signature, and the elon- 
gated correlation functions corresponding to the other two 
cluster samples in this figure. This elongation is most se- 
vere in the case of Abell clusters, as it has already been 
pointed out in several previous works (Bahcall & Soneira 
1983; Postman, Huchra & Geller 1992), and has been rec- 
onciled by considering the effects of cluster alignments in 
supercluster structures (Bahcall, Soneira & Burgett 1986), 
and the likely inclusion of spurious clusters in the sample 
due to projection effects (Sutherland 1988; Sutherland & 
Efstathiou 1991; see also Lucey 1983). APM clusters show 
a much milder elongation, which can be easily reproduced 



in numerical simulations when a rodsliift error bigger than 
700km/s is added to the mock cluster positions (Padilla & 
Baugh 2002) although contamination can not be totally ex- 
cluded. 

In an attempt to explore and quantify the origin of the 
problem of large elongations of the clustering signal we study 
the effects of errors in cluster distances by considering dif- 
ferent number of redshifts used in the determination of clus- 
ter distances, as well as the possible projection effects. In 
order to do this, wo use results obtained from mock cata- 
logues where clusters of galEixies are identified from angular 
data (see Paper I). Such a comparison will also be useful 
in providing a quantitative determination of the degree of 
projection effects present in the Abell cluster sample. 

3.1.1 Effects of different number of redshifts 

In order to assess to which extent the use of samples of clus- 
ters with large number of rcdsliift measurements, Uz, reduces 
the distortion of the clustering pattern, we have explored 
samples with different values of n^. 

Figure 2 shows the results from clusters whose distances 
have been obtained using more than ten galaxy redshifts, 
> 10 (left panel), and results for clusters with Uz < 10 
(right hand side panel) . The net effect of having a large riz is 
to lower the amplitude of the line of sight elongation. It can 
be seen that, even though this effect is not small, it is insuf- 
ficient to reproduce the fiattening observed in less massive 
groups of galaxies as seen in figure 1. The results shown in 
the two panels in figure 2 can be compared to figure 6, paper 
I, for different riz values. These samples were constructed in 
order to match as closely as possible the procedure followed 
to construct cluster samples obtained from angular informa- 
tion. By comparison, the trend of increasing anisotropics for 
lower number of redshift measurements used to determine 
cluster distances is present in both, the observational and 
simulated samples. 

An interesting point to be noted is the clear enhance- 
ment of the correlation length in the a direction for the small 
TT bin (which we will call ao from now on) in the observa- 
tional and mock samples with low , with respect to those 
with more redshift measurements. By inspection to figure 
8 we see that this difference is also present in the redshift- 
space correlation length sq. This is in agreement with more 
general results obtained from mock cluster catalogues iden- 
tified from angular simulated data. This can be seen in figure 
10 of Paper I, where a decrement in so for higher values of 
Uz can be easily seen. This indicates that it is likely that the 
difference in correlation amplitudes indicated by the value 
of So shows the diminished effects of projection biases in the 
high Hz Abell sample, rather than a difference in average 
cluster mass. 

We will come back to this point in the following sections, 
where we will examine subsamples of optically and x-ray 
selected clusters of galaxies. 

3.1.2 Results from different from Bautz-Morgan type 
samples 

We have also analysed samples of Abell clusters with differ- 
ent Bautz-Morgan (BM) types. The BM type is indicative 
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Figure 1. Correlation function of UZC groups (left hand side panel), APM clusters (Dalton et al. 1992; middle panel), and the full 
sample of Abell clusters with richness class r > 1 and measured redshifts (Struble & Rood, 1999, right panel). The transitions between 
different shadings correspond to fixed values of ^ = 0.6, 0.8, 1.0, 1.2, and 1.4 levels, which are the values used later to infer the relative 
velocities, ■wi2. The thick line corresponds to the ^ = 1 level, and the thin lines show the expected contours for a spherically symmetric 

correlation function. 




Figure 2. Correlation function of samples of Abcll with different number of redshifts used in their distance determinations. The left 
panel shows correlation function contours obtained from Abell clusters whose distances have been obtained using more than ten galaxy 
redshifts. The results for Abell clusters with less than ten redshift measurements are shown on the right hand side panel. Shadings and 
line conventions are as in figure 1. 



of the regularity of the distribution of galaxies in a cluster; 
the larger the BM type, the more irregular the shape of the 
cluster is. 

The results from studying samples of different BM type 
can be seen in figure 3, where the left panel shows the corre- 
lation contours for clusters with Bautz-Morgan type I, I-II, 
and II, corresponding to more regular clusters, and the right 
panel shows results from BM types II-III and III, that is the 
most irregular ones. The particular choice of BM type to 
define the samples corresponds to a rough classification of 
regular and irregular clusters. Even though the differences 
between both panels axe not conclusive, there is a larger 
elongation in the correlation function of clusters with larger 



BM types (the inferred relative velocities for these cluster 
samples can be seen in figure 8). 

The larger elongation seen in the correlation function 
of high BM type clusters in figure 3 indicates that they are 
slightly more prone to the inclusion of spurious background 
or foreground galaxies than low BM types, or even to be 
results of projection effects in the angular distribution of 
galaxies. These problems seem also to be present in the low 
BM type clusters, but to a lesser extent. 

It can also be appreciated that the estimated correlation 
length (To has similar values for the low and high BM type 
clusters (figure 3). This fact would suggest a similar real- 
space correlation length for both samples. However, we ex- 
pect that the most massive clusters would show a more reg- 
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Figure 3. Rcdsliift-spacc correlation function for Abell Clusters with different Bautz-Morgan type. The left panel shows the correlation 
contours for clusters with low Bautz-Morgan type (I, I-II, and II). Results from high Bautz-Morgan types (II-III and III) are shown in 
the right hand side panel. Shadings and line conventions are as in figure 1. 



ular galaxy distribution, and therefore, be of low BM types, 
whereas high BM typo clusters would correspond to less 
massive clusters, with lower galaxy richness counts, more 
irregular shapes, and lower correlation lengths. Since high 
BM clusters axe likely to be subject to larger projection ef- 
fects due to their irregular appearance they could have a 
biased large estimate of ao as shown in paper I, clarifying 
the results of ao seen in this figure. 

We next measure the redshift-space correlation length. 
So, for both sub-samples. In this case, we obtain a higher 
correlation length for the low BM type clusters (see figure 
8), in support of the presence of a correlation between mass 
and BM type proposed in the last paragraph. 

Furthermore, as we discussed in the previous sub- 
section, it can be argued that both ao and so depend on 
the degree of projection effects present in a sample. The 
analyses carried out in Paper I indicate that the effects of 
increasing Uz on so, for samples of clusters identified from 
angular data, is that of making so ever smaller and closer 
to its actual value. This result implies that the fact that we 
obtain a larger value of so for clusters with low BM type 
actually indicates that this sample contains more massive 
clusters than that with high BM type clusters. 

Prom these results, we notice that we can reduce the 
degree of contamination in a sample of Abell clusters by 
considering subsamples of regularly shaped systems. How- 
ever, it should be noticed that projection effects cannot be 
completely removed in the Abell catalogue by considering 
subsamples with low BM type nor large Uz . 

3.2 X-ray selected clusters of galaxies 

In order to provide an insight to the characteristics of cluster 
samples with X-ray selection criteria, we study two samples 
of X-ray selected clusters of galaxies, the X-ray Bright Abell 
Cluster Sample (XBACS, Ebeling et al. 1996) and the Ex- 



tended ROSAT Brightest Cluster Sample (EBCS, Ebeling 
et al. 2000). The fact that these clusters have been identi- 
fied from the X-ray emission of hot gas trapped in the huge 
potential wells of massive clusters of galaxies, makes these 
samples more reliable than cluster samples selected from 
2D data. This is a consequence of the fact that the emission 
from the hot gas is proportional to the square of the den- 
sity, which makes it a lot more difficult to have a spurious 
cluster of galaxies resulting from the alignment of small hot 
gas clouds along the line of sight. 

Clusters of galaxies in the XBACS sample, which com- 
prises Abell clusters whose emission in X-rays has been con- 
firmed, would be less affected by projection effects present in 
the Abell sample, since only real clusters would show hot gas 
emission. On the other hand, in the EBCS sample, sources 
were identified from an X-ray survey and then galaxy red- 
shift distributions were analysed to find the optical counter- 
parts. These clusters should be largely free from projection 
effects, and should not include identification bieises associ- 
ated to the Abell cluster finding algorithm. However, the 
reliability of X-ray cluster distances obtained from correla- 
tions with redshift surveys are subject to uncertainties that 
could affect the results. The correlation functions measured 
for the EBCS and XBACS samples are shown in the left 
and right hand side panels of figure 4 respectively. As can 
be seen, the EBCS sample shows only a slight elongation 
along the line of sight which can be easily reconciled with 
the presence of distance measurement errors similar to those 
necessary to explain the anisotropics observed in the APM 
correlation function (Padilla & Baugh, 2002). On the other 
hand, the line of sight distortion of the XBACS correlation 
function shown in the right panel of this figure provides a 
strong indication of the presence of projection effects. Given 
their larger correlation amplitude, XBACS should be subject 
to a strong infall signal, which is not observed. The results 
from Paper I indicate that the projection effects seen in the 
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Figure 4. Correlation function contours from EBCS (left hand side panel) and XBACS clusters (right hand side panel). Shadings and 
line conventions are as in figure 1. 



correlation function of mock clusters identified from angu- 
lar data (sample 3) are still present in a set of clusters with 
angular positions coincident with clusters identified in 3D 
(sample 4). The studies performed in Paper I show that the 
elongations in ^ (cr, tt) found in this mock cluster catalogue 
are quite severe, and are mainly produced by the erroneous 
use of foreground and background galaxies in assigning the 
distances to clusters. 



3.2.1 Luminosity effects in x-ray selected samples 

Even though the elongation of the correlation function con- 
tours along the line of sight present in the total EBCS sam- 
ple can be accounted for by redshift measurement errors, 
wo explore EBCS subsamplcs defined by limiting the range 
of intracluster gas luminosity in x-rays. We define a limit 
Liimit = 3 X 10^^ VK which divides the total sample of EBCS 
clusters into two subsamples with similar number of clusters. 
The resulting iso-correlation levels can be seen in figure 5, 
where the left panel shows the results for Lx < Li,m?t; and 
the right hand side panel, those for > Lumit- This is a 
somewhat unexpected result, since high luminosity clusters 
show a very elongated pattern, whereas the least luminous 
clusters show the expected flattened pattern indicative of 
less significant projection effects. 

A possible source for this effect could rely on contamina- 
tion by foreground structures which are likely to affect more 
strongly clusters at greater distances. Those contamination 
effects could bias distance estimates and the X-ray luminos- 
ity. In figure 6 we show the distribution of cluster distances 
for the high and the low x-ray luminosity subsamples which 
support this scenario. 

The resulting values of relative velocities, rcdsliift-spacc 
correlation lengths and bias factors corresponding to the 
EBCS subsamples defined by the luminosity cut Lumu are 
presented in figures 8 and 9 respectively. As can be seen, 
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Figure 6. Redshift distribution for the two EBCS subsamplcs 
defined by a luminosity cut Lx < iiimit (upward hatched region), 
and iamzt < (downward hatched region). 



low x-ray luminosity EBCS clusters present a remarkably 
low pairwise velocity dispersion, comparable to that result- 
ing from the analysis of UZC groups, identified from redshift 
data. 



3.2.2 Subsamples of different BM type effects in x-ray 

clusters 

In this subsection we analyse the redshift space clustering 
pattern arising in samples of x-ray clusters with different 
BM types. This is aimed to study the combination of projec- 
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Figure 5. Correlation function for subsamples of EBCS clusters with different luminosities. The left panel shows the correlation contours 
for clusters with Lx < 3x lO^'' W. Results from using Lx > 3x 10^'^ W axe shown in the right hand side panel. Shadings and line conventions 
are as in figure 1. 
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Figure 7. Correlation function for subsamples of XBACS clusters with different Bautz-Morgan type. The left panel shows the correlation 
contours for clusters with Bautz-Morgan type I, I-II and II. Results from higher Bautz-Morgan types are shown in the right hand side 
panel. Shadings and line conventions are as in figure 1. 



tion effects which could manifest in the correlation function. 
Since BM types are available for Abell clusters, we restrict 
the analysis to the XBACS sample. 

The results can be seen in figure 7, where tfic left 
panel shows tfie correlation contours for clusters with Bautz- 
Morgan type I, I-II, and II and the right panel shows results 
from BM types II-III and III. The choice of BM types for the 
different panels is the same as that used for Abell clusters, 
but the difference in elongations is more readily apparent. 
The figure shows a significantly larger elongation for the 



sample of clusters with larger BM types, as was the case 
with Abell clusters of different BM type. The larger corre- 
lation length (To for the low BM type clusters, which can 
also be seen as a larger value of so (see figure 8) can be 
appreciated once more. This supports the earlier hypothe- 
sis that the most massive clusters would be more frequently 
identified as low BM type clusters, and therefore show a 
large correlation length. This figure also supports the idea 
that the larger elongations in the high BM type clusters 
are probably produced by a larger number of background or 
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Figure 9. Bias factors for the observational cluster samples stud- 
ied in this paper. The circles with errorbars show the results from 
minimising equation 2, and the solid line shows the predicted ef- 
fective bias from the CDM model and the Sheth, Mo & Tormen 
(2001) mass function. The range of a<;ceptable values for the ef- 
fective bias factors axe shown as dotted lines. The ordering of the 
samples is the same as in figure 8. 



Figure 8. Pairwise velocities (upper panel) and correlation 
lengths (lower panel) obtained for the different observational sam- 
ples. The samples are distributed along the x-axis, such that their 
values of wi2 decrease to the right of the plot. This ordering is 
maintained in the lower panel. 

foreground galaxies used in assigning distances to high BM 
type clusters. 

3.3 Relative velocities and correlation lengths in 
observational cluster samples 

In order to summarise the results obtained from the study of 

observational samples of clusters of galaxies identified from 
angular data, we present the values of relative velocities and 
redshift-space correlation lengths obtained from the correla- 
tion functions shown in previous figures. 

We find the relative velocities by minimising equation 
2. The real space correlation function used in this equa- 
tion corresponds to a CDM power spectrum with parameters 
^matter = 0.3, Qa = 0.7, a CDM shape parameter F = 0.2, 
and (78 = 0.9. 

Based on the results of the study of galaxies and groups 
of galaxies, and supported by simulations of the hierarchi- 
cal processes thought to drive the formation of structure, 
we expect to find smaller relative velocities for clusters of 
galaxies than are found for galaxies and groups. The upper 
panel of figure 8 shows relative velocity on the y-axis, and 
observed sample on the x-axis. The errors in relative ve- 
locity result from the different estimates obtained from dif- 
ferent correlation function levels, and show 1-cr confidence 
levels. The ordering of the samples along the x-axis is such 



that relative velocities decrease to the right. If the cluster 
samples wore free from projection effects of any kind, this 
ordering of samples would result in showing groups on the 
left and the most massive clusters on the right. However, 
as can be seen, every cluster sample, with the exemption of 
EBCS clusters with low x-ray luminosity, is located to the 
left of this plot, showing their large contamination due to 
projection effects. The lower panel shows the rcdshift-space 
correlation length with errorbars obtained from performing 
a power law fit to the measured redshift-space correlation 
functions over a range of scales dependent on the sample 
analysed Tmin < r < rmax, where 2.5 < rmin/h~^Mpc< 10, 
and 28 < rmax /h~^Mpc< 80 (the ordering of samples in this 
panel is the same as in the previous one). There are signs 
of a weak correlation between the amplitude of the relative 
velocities and the correlation length for the cluster samples. 
According to the results from Paper I, the degree of pro- 
jection effects correlates with both, the value of the relative 
velocities and the correlation length, so a similar trend in 
both parameters was expected; the deviations from this be- 
haviour reflect the different minimum mass and X-ray flux 
thresholds in the different observational samples, which in 
the case of the mock samples with different Uz was always 
the same by construction. 

In a hierarchical scenario of structure formation it is 
expected that cluster relative velocities should not be larger 
than galaxy pairwise velocities (~ 4:00km/ s). The observa- 
tions suggest apparent large pairwise velocities W12, an effect 
which could have a large contribution from redshift measure- 
ment errors, and even more important contributions from 
large systematics due to projection biases present in some 
of the samples identified from angular data. 
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This is the case with Abell clusters and the XBACS 
sample, except when considering subsamples of clusters with 
low BM type. Large relative velocities are inferred for lumi- 
nous EBCS clusters, giving hints for the presence of con- 
tamination duo to projection effects. However, the results 
for EBCS of low x-ray luminosity, show significantly smaller 
pairwise velocities comparable to those found for the UZC 
groups, and consistent with the predictions of a hierarchical 
clustering scenario. 

3.4 Bias factors 

The comparison between the measured and predicted 
redshift-space correlation function can be used to obtain 
estimates of the bias parameter of the cluster sample un- 
der consideration, using equation 2. In figure 9 we show the 
different bias parameters obtained for the different samples 
analysed. The errors on the bias parameter come from the 
different results obtained from the comparison of 5 levels of 
correlation function. 

These values can be compared to what is oxpoctod from 
a particular CDM model, using the expression for an effec- 
tive bias as presented in Padilla & Baugh (2002), which is 
a weighted average of the bias of a sample of clusters of 
mass M, derived by Sheth, Mo & Tormen (2001), and re- 
quires knowledge of the space density of the cluster sample. 
The results from using this equation are valid if the cluster 
sample is complete above a minimum mass threshold. 

We have derived approximate estimates of the num- 
ber density of clusters by assuming that the cluster sam- 
ples arc complete above a minimum cluster mass limit. This 
rough approximation only applies to optically selected clus- 
ters since x-ray selected samples are flux-limited with a min- 
imum mass which strongly depends on redshift. Then we cal- 
culate the space density of optical cluster samples by simply 
counting the number of objects up to a distance where the 
redshift distribution departs from the behaviour due to 
the influence of a selection function. The effective CDM bias 
for the measured cluster number densities are then derived 
straightforward following Padilla & Baugh (2002). Given the 
uncertainties in the assumptions mentioned before, and the 
large scatter between cluster richness and virial mass, we 
just provide a suitable upper limit for the bias parameter, 
namely the largest value obtained, b = 3.5, corresponding 
to the Abell sample with low BM type (this is equivalent to 
a mean intercluster separation dc — 60h~^Mpc or a number 
density n ~ 5 x 10~®h^Mpc~^). This upper limit is shown 
as a dotted line in figure 9. 

For some of the samples studied, the resulting bias pa- 
rameters obtained from the minimisation are well above 
our estimated upper limit (figure 9). However, UZC groups, 
APM clusters, EBCS total sample and low x-ray luminosity 
subsample, and Abell clusters with high riz , show lower val- 
ues of the bias parameter, compatible with our expectations. 

4 AN INSIGHT ON THE RESULTS FROM 
GROUPS AND CLUSTERS IDENTIFIED 
FROM REDSHIFT DATA 

The UZC group sample was constructed by identifying 
galaxy systems from the 3D distribution of galaxies in red- 



shift space from the Updated Zwicky catalogue. By inspec- 
tion of figure 8, the great leap between the results of UZC 
groups relative velocities, and those obtained from samples 
of clusters identified from angular data, cither in the optical 
or x-ray wavelength range, is clear. This is also apparent in 
figure 1, where the anisotropies present in the correlation 
function of the UZC groups show markedly flattened iso- 
correlation contours, a pattern we almost completely failed 
to reproduce by measuring the correlation function of other 
cluster samples throughout this paper with the exception of 
the low x-ray luminosity EBCS subsample. 

In this section, we explore in more detail the relative ve- 
locities, redshift-space correlation lengths and bias parame- 
ters obtained from the correlation function of UZC groups 
as a function of group richness counts. We remark the differ- 
ence between the number of member galaxies in a group and 
the values of riz from the previous analysis, for in this case 
the number of galaxies is associated to the identification pro- 
cess and not only to the number of redshift measurements 
in the field of a cluster. 

We note that the number of member galaxies is cor- 
related to some extent with the mass of the groups. This 
can be seen in figure 10, where we show the distribution of 
group virial masses for two different upper limits in number 
of member galaxies. The normalisation constant by which 
both distribution functions were multiplied, is set up so as 
to make the integral of the mass distribution function of 
groups with rimember < 5 cqual to 1. It can be seen that 
the sample with a low upper limit for rimember (dark shaded 
area) shows a marked relative lack of high mass groups when 
compared to the sample with a high upper limit in member 
galaxies (light shaded area). One of the effects this leads 
to is the intrinsic change in the clustering amplitude of the 
different rimember subsamples. It should be noted that sam- 
ples with different riz in mock sample 2, are not affected by 
this bias. This is readily apparent when recalling a result 
from Paper I, where it was shown that this sample is con- 
sistent with a constant redshift-space correlation length for 
any value of . 

Following the procedures given in sections 3.3 and 3.4, 
wo calculate the relative velocities, redshift-space correlation 
functions, and bias corresponding to subsamples of the UZC 
groups with different number of member galaxies. The re- 
sults can be seen in figure 11. The upper panel shows values 
for the relative velocities which are consistent with the pub- 
lished value of wi2 = 250±110km/s (Padilla et al. 2001) for 
any value of rimember • There are also signs of a somewhat 
larger value of wi2 for small rimember , but within the size of 
the crrorbars. The middle panel shows the values of so as 
a function of rimember- As can be seen, there is a tendency 
for higher values of so as the number of member galaxies in- 
creases. This increase is also present in the measurements of 
bieis factors, shown in the lower panel, and is consistent with 
the results from figure 10, which show an increasingly larger 
average mass for samples with larger rimember- The higher 
the average mass of a sample, the higher its correlation am- 
plitude and bias factor will be. The increase in so and bias 
factor are only noticeable for rimember < 10. Most of the 
samples with a larger number of member galaxies show sta- 
ble and constant results for the three parameters in study, 
wi2, So and bias factor, in good agreement with the results 
obtained from mock sample 2, which also showed this be- 



© 0000 RAS, MNRAS 000, 000-000 



-1 - 



o 

(U 

u 

"2 -2 

so 





T — I — I — I — I — I — I — I — I — I — 1 — 1 — 1 — 1 — 1 — r 



Anisotropics in the redshift- space correlations II: Observations 

500 
400 



11 



' ' ' I ' ' ' 
UZC Groups 



nmember<40 



^ member ''^^ 



12 



13 14 
log,o(M/mg) 



15 



Figure 10. Virial masses distribution function for the UZC 
groups, for different subsamples, defined by an upper limit in 
i^member- The light shaded area corresponds to TimemiieT- < 40, 
and the dark shaded area to n^nember < 5- 



haviour. The errorbars in figure 11 shown for W12 and bias 

factor are obtained from the different results obtained by us- 
ing five correlation function amplitudes = 0.6, 0.8, 1.0, 1.2, 
and 1.4). 

The main conclusion of this section is the remarkable 
stability of the results for the UZC groups. This provides a 
clear indication that redshift space identification of clusters 
is a powerful procedure to obtain samples with no signif- 
icant contamination. The weak dependence of the results 
on a wide range of number of member galaxies reflects the 
reliability of the identification procedure. 



5 CONCLUSIONS 

We have explored the correlation function of clusters in red- 
shift space for different observational samples of galaxy clus- 
ters. We have studied the consequences of different obser- 
vational issues on the distortion of the clustering pattern, 

such as different sources of cluster identification, number of 
redshift measurements, cluster morphology and X-ray lumi- 
nosity. 

Padilla et al. (2001) showed that the results from groups 
of galaxies identified from 3-dimensional information avail- 
able in the UZC, arc consistent with the hierarchical cluster- 
ing paradigm which predicts that we should be able to find 
objects that arc still falling towards more massive struc- 
tures. This can be observed in the anisotropy of the correla- 
tion function contours of the UZC groups, which shows the 
fiattening produced by this infall motion. This is not easily 
detectable for galaxies since they are in a strongly non-linear 
regime and the pairwise peculiar velocities dominate the dy- 
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Figure 11. UZC groups relative velocities (upper panel), 
redshift-space correlation lengths (middle panel), and bias fac- 
tors (lower panel) as a function of minimum number of member 
galaxies. 



The automatic cluster finding algorithms used in the 
identification of APM clusters, make them a more reliable 
sample of clusters, which have been shown to be largely 
free of projection effects. Although contamination can not 
be totally excluded, the inclusion of an error in the APM 
cluster distance measurement is sufficient to account for for 
the elongations found in the APM clusters redshift-space 
correlation function (Padilla & Baugh 2002). 

On the other hand, Abell clusters show clear signs of 
projection effects (this work, Sutherland 1988; Sutherland 
& Efstathiou 1991; see also Lucey 1983), which can be re- 
duced by restricting the sample of clusters in richness class 
(Miller et al. 1999), for example. In this work, we found that 
the elongations in the correlation function isopleths can be 
reduced when restricting the sample to have a low BM type, 
which is a measure of the regularity of the galaxy distribu- 
tion in a cluster. Wc have also shown that this could be 
related to an increase in the richness of the clusters consid- 
ered, since our correlation function measurements show an 
increase in the correlation length. This makes sense, given 
that it is reasonable to suppose that richer clusters of galax- 
ies have a more regular shape. In the end, the anisotropics 
left after the restriction of the Abell sample cannot be easily 
reconciled with errors in distance measurement from using a 
small number of galaxies. The results indicate the presence 
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of galaxies not physically bounded to the cluster, and also a 
fraction of spurious clusters resulting from structures along 
the line of sight. 

This situation is not significantly improved in the case of 
X-ray clusters. Even though the XBACS sample is composed 
of Abell clusters with confirmed emission in X-rays, they 
also show significant elongations along the line of sight when 
compared to the results from the Abell sample. We found 
that restricting the sample of XBACS clusters to have a 
low BM type produces a decrease in the elongation of the 
iso-correlation contours along the line of sight. 

On the other hand, EBCS clusters, which are iden- 
tified primarily as X-ray sources, show an interesting be- 
haviour that deserves further attention. Low X-ray luminos- 
ity clusters, L^; < 3 x show little correlation function 
anisotropy indicating the lack of significant projection ef- 
fects. Clusters of high X-ray luminosity, on the other hand, 
exhibit strong line of sight distortions, indicating the pres- 
ence of large biases in this subsample. It can be argued 
that luminous clusters are more distant and therefore most 
strongly affected by contamination effects that could bias 
distance and x-ray luminosity estimates. 

The results shown in figures 8 and 9, which correspond 
to observational datasets identified from angular data, in- 
dicate that the values of relative velocities obtained from 
cluster samples are usually well over the values preferred by 
hierarchical clustering. It is also noticeable the large varia- 
tions in redshift-space correlation lengths which, as shown 
in Paper I, is also infiuenced by the level of contamination 
present in the sample. Those observational samples less af- 
fected by projection effects, show values of the bias paxam- 
eter compatible with a suitable upper limit inferred using a 
number density threshold n — 5x 10^''h'^Mpc^'' (equivalent 
to a mean intercluster separation dc ^ 60h~^Mpc), and a 
ACDM model in combination with the Sheth, Mo & Tormen 
(1999) mass function. 

With the advent of surveys of the next generation, the 
construction of new group and cluster catalogues from deep 
spectroscopic data based on accurate photometry, will pro- 
vide more reliable and well controlled cluster samples. In 
the light of our results for presently available samples, we 
find that detailed analysis of the redshift-space correlation 
function anisotropies can give a deep insight on sample char- 
acteristics and invaluable information on their implications. 

ACKNOWLEDGMENTS 

This work was supported in part by CONICET, Argentina, 
and the PPARC rolling grant at the University of Durham. 
DGL acknowledges support from the John Simon Guggen- 
heim Memorial Foundation. We thank the Referee for invalu- 
able comments and advice which greatly improved the previ- 
ous version of the paper. We have benefited from helpful dis- 
cussions with Carlton Baugh. We acknowledge the Durham 
Extragalactic Astronomy Group and the Virgo Consortium 
for making the Hubble Volume simulation mock catalogues 
available. 



Bahcall, N.A., Cen, R.Y., Gramann, M., 1994, ApJ, 430, L13. 
Bahcall, N.A., Soneira, R.M., 1983, ApJ, 270, 20. 
Bahcall, N.A., Soneira, R.M., Burgett, W.S., 1986, ApJ, 311, 15. 
Bautz, L.P. & Morgan, W.W., 1970, ApJ, 162, 149. 
Borgani, S., Plionis, M., Kolokotronis, V., 1999, MNRAS, 305, 
866. 

Cole, S., Hatton, S.J., Weinberg, D.H., Prenk, C.S., 1998, MN- 
RAS, 300, 945. 

Collins, C.A., Guzzo, L., Bochringcr, H., Schucckcr, P, Chincar- 
ini, G., Cruddace, R., De Grandi, S, MacGillivray, H.T., Neu- 
mann, D.M., Schindler, S., Shaver, P., Voges, W., 2000, MN- 
RAS, 319, 939. 

Croft, R.A.C., Efstathiou, G., 1994, MNRAS, 268, L23. 

Croft, R.A.C., Dalton, G.B., Efstathiou, G., Sutherland, W.J., 
Maddox, S.J., 1997, MNRAS, 291, 305. 

Dalton, G.B., Efstathiou, G., Maddox, S.J., Sutherland, W.J., 
1992, ApJ, 390, LI. 

Dalton, G.B., Efstathiou, G., Maddox, S.J., Sutherland, W.J., 
1994, MNRAS, 269, 151. 

Dalton, G.B., Maddox, S.J., Sutherland, W.J., Efstathiou, G., 
1997, MNRAS, 289, 263. 

Davis, M., & Peebles, P.J.E., 1983, ApJ, 267, 465. 

Ebeling, H.; Voges, W.; Bohringer, H.; Edge, A. C; Huchra, J. 
P.; Briel, U. G., 1996, MNRAS, 281, 799. 

Ebeling, H.; Edge, A. C.; Allen, S. W.; Crawford, C. S.; Fabian, 
A. C; Huchra, J. P., 2000, MNRAS, 318, 333. 

Eiscnstein, D., & Hu, W., 1998, ApJ, 496, 605. 

Evrard, A. E.; MacFarland, T. J.; Couchman, H. M. P.; Colberg, 
J. M.; Yoshida, N.: White, S. D. M.; Jenkins, A.; Frenk, C. 
S.; Pearce, F. R.; Peacock, J. A.; Thomas, P. A. 2002, ApJ, 
573, 7. 

Huchra, J.P., & Geller, M.J., 1982, ApJ, 256, 423. 
Lucey, J. R., 1983, MNRAS, 204, 33. 

Lumsden, S.L., Nichol, R.C., Collins, C.A., Guzzo, L., 1992, MN- 
RAS, 258, 1. 

Miller, C.J., Krughoff, K.S., Batuski, D.J. & Hill, J.M., 2002, AJ, 
124, 1918. 

Miller, C.J., Batuski, D.J., Slingend, K.A., HiU, J.M., 1999, ApJ, 

523, 492. 

Netterfield, C. B. et al., 2002, ApJ, 571, 604. 
PadiUa, N.D., & Lambas, D.G., 2003, MNRAS, submitted (Paper 
I)- 

Padilla, N.D., & Baugh, CM., 2002, MNRAS, 329, 431. 
Padilla, N.D., Merchan, M.E., Valotto, C.A., Lambas, D.G., 

Maia, M.A.G., 2001, ApJ, 551, 873. 
Peacock, J. A. & West, M.J., 1992, MNRAS, 259, 494. 
Postman, M., Huchra, J. P., Geller, M.J., 1992, ApJ, 384, 404. 
Sheth, R.K., Mo, H.J., Tormen, G., 2001, MNRAS, 323, 1. 
Struble, M.F. & Rood, H.J., 1999, ApJS, 125, 35. 
Sutherland, W., 1988, MNRAS, 234, 159. 
Sutherland, W.J., Efstathiou, G., 1991, MNRAS, 248, 159. 
van Haarlem, M.P., Prenk, C.S., White, S.D.M., MNRAS, 1997, 

287, 817. 



REFERENCES 

Abadi, M.G., Lambas, D.G., Muriel, H., 1998, ApJ, 507, 526. 



© 0000 RAS, MNRAS 000, 000-000 



